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both the minimal and the non-minimal SUGRA models can generate neutralino- 
proton cross-sections at the level compatible with the signals reported in the DAMA 
experiment at the Gran Sasso National Laboratory. Effects of proton stability on the 



(— I \ analysis of the DAMA data in the minimal and the non-minimal SUGRA models are 



also discussed. 



1. INTRODUCTION 



Last year, the DAMA experiment examined the possibility of the direct detection of 
Milky Way wimps using the annual modulation signal. Based on a set up of approximately 
100 kg of radiopure NaI{Ti) detectors in the Gran Sasso National Laboratory and 4549 kg- 
day of data, they found an indication of such a signal. More recently, an additional 14,962 
kg-day of data has significantly strengthened the statistical analysis favoring the presence 
of a yearly modulation signalH with a wimp mass and proton cross section of 

M^ = {59tll)GeV; ^a^.p = (7.0^?;^) x 10"^ (1) 



where ^ = p«,/po, Pw is the local Milky Way wimp mass density and po = O.SGeV^cm"^. 

In the following, we will analyze the possible consequences of such a signal within the 
framework of supergravity grand unification models with R-parity invariance and gravity 
mediated supersymmetry (SUSY) breakingH'Q. Such models automatically predict the ex- 
istance of cold dark matter (CDM) in the universe, i.e. the relic lightest supersymmetric 
particle (LSP) remaining from the Big Bang. Further, over most of the SUSY parameter 
space, the LSP is the lightest neutralino, Xi; and the calculated relic density of Xi is in 
accord with astronomical estimates over a significant part of the parameter space. For the 
relic density, we will assume the range 0.05 < Vt^ah"^ < 0.30, where Vt^o = p^o/pc, where p^o 
is the mean relic density in the universe, pc = SH^/SttGn, G^ = Newton's constant, and 
the Hubble constant H is parameterized hj H = {100kmsec~^Mpc~^)h. 

The annual modulation effect arises due to the motion of the Earth around the Sun. 
Thus, ve, the velocity of the Earth relative to the Galaxy is ve = vs + vq cos 7 cos (^^(t — to) 
where vs is the Sun's velocity relative to the Galaxy {vs = 232km/ s), Vq is the Earth's 
orbital velocity around the Sun {vq = 30km/ s) and 7 is the angle of inclination of the plane 
of the Earth's orbit relative to the galactic plane (7 = 60°). One has u = 2tt/T (T = 1 year) 
and the maximum velocity occurs at day to = 155.2 (June 2). The change in the Earth's 
velocity relative to the incident wimps leads to a yearly modulation of the scattering event 
rates of about 7%. 

The calculation of the xi — p cross section proceeds as foUowsS. One first calculates the 
relic density of Xi^ limiting the SUSY parameter space so that the above constraints on 
Q^oh'^ are obeyed. Within this constrained parameter space, we then calculate the Xi ~ P 
cross section for incident halo Xi on the terrestial target. In comparing the theoretical cr^o_p 
with the data, a number of uncertainties arise due mainly to the lack of knowledge of input 
parameters. We estimate an error of a factor of 2-3 in these calculations. In addition, p^o 
may vary from about (0.2 — 0.7)GeVcm^'^ (i.e., 0.7 ~ ^ ~ 2.3) . 

Supergravity models have a wide range of applicability. Thus, once one phenomena 
begins to fix the SUSY parameters, it effects predictions in other areas. We will examine 



here the effects the DAMA data has on SUSY mass spectrum predictions at future accelerator 
searches. In addition, most supergravity models predict the existence of proton decay. While 
the predictions for proton decay are more model dependent than other parts of the theory, 
there is a strong correlation in the SUSY parameter space between DM detector event rates 
and the expected proton lifetimeB, and we will discuss this below. 

2. SUPERGRAVITY MODELS 

We consider supergravity (SUGRA) models where supersymmetry is broken in a hid- 
den sector at the Planck scale [Mp = 2.4 x lO^^Gel^) by supergravity interactions and 
transmitted to the physical sector by supergravity, giving rise to soft breaking termsB. If 
the hidden sector interactions are generation independent, one obtains the simplest model, 
mSUGRA, with universal soft breaking parameters at the GUT scale Mq ~ 2 x lO^^GeV. 
The renormalization group equations (RGE) then show that the SUSY soft breaking at Mq 
generates SU{2) x f/(l) breaking at the electroweak scale. This model then depends on four 
parameters and one sign (in addition to the parameters of the SM) . One can choose these 
parameters to be the following: rrio, the universal scalar soft breaking mass at Mq] fni/2, 
the universal gaugino mass at Mq (or alternately the gluino mass nig = {a^{Mz) / aG)'mi/2 
where ac — 1/24 is the coupling constant at Mc); Aq, the universal cubic soft break- 
ing parameter (or alternately At, the t-quark parameter at the electroweak scale); and 
tan/3 = < H2 >/< Hi > where < H2 > gives rise to the up quark masses and < iJi > to 
the down quark and the lepton masses. The RGE then determines fi^ (where /i is the Higgs 
mixing parameter in the superpotential term fiHiH2), leaving the sign of /i arbitrary. 

The mSUGRA model contains relatively few new parameters, and existing data has 
begun to limit this parameter space. Thus, the measured value of the t-quark mass and the 
b —>■ s + '-f branching ratio eliminate most of the /i < and At < part of the parameter 
space0. However, the supergravity formalism allows for non-universal soft breaking^. While 
the universality of the soft breaking mass mo guarantees the suppression of flavor changing 



neutral currents (FCNC), the FCNC are not sensitive to the Higgs mass or to the third 
generation non-universahties and string models can allow for such non-universality. We 
parameterize this possibility by the following Higgs soft breaking masses at M^, 

mjj, = m°(l + 5i); ml^ = ml{l + 62) (2) 

and the following third generation sfermion masses: 

ml^ = ml{l + 63)] ml^ = ml{l + 64); ml^ = ml{l + 65) 



m 



L = ^0(1 + Se); ml = ml{l + 6r) (3) 



where qi = {tL,bL), ur = tR,iL = {pl.^l), etc. uiq is the universal soft breaking mass of 
the first two generations, and 5i represents the non-universal deviations. For soft breaking 
occurring above M^, and for GUT groups having an SU(5) subgroup (e.g. SU(N), A^ > 5, 
SO(N), A^ > 10, Eq, etc.) with matter embedded into 10 and 5 representations in the usual 
way, one has ^3 = 54 = ^5 and 5^ = 67. We assume in the following that | 5j |< 1. 

The non-universal corrections enter sensitively in fi^. For tan/5 < 25,^5,56 and 67 make 
only small contributions, and a closed form expression can be obtained for /z^ at the elec- 
troweak scaled 



2 t ..1 — 3Dq 1 1 — Dq 1 + Dq 1 e m 2 

^ = ^^3y[( 2 ^2) + (^ — ("^3 + ^4) — S2 + ^5i)]mo 



+ ^T^[-(l - I?o)— + CX/o] - -Ml + —^^5,(1 - ^i^) + 1 loop terms (4) 

where t = tan/3, C^, = iDo(l - Do){H3/Fy + e - g/t^, and Dq ~ 1 - m2/(200sin/3)2. 
Dq vanishes at the t-quark Landau pole (for rrit = 175GeV, Dq < 0.23) and Ar = At — 
m.i/2{H2 — H3/F) is the residue at the Landau pole {Ar = At — 0.61(03/00)^-1/2), i.e. 
Aq = Ar/ Dq — {H3/ F)mi/2- Sq = TrYm? where Y is the hypercharge and m? are the 
squark and Higgs masses at Mq of Eqs. (2,3). The form factors H2-, H^, F, e, g are given in 
Ref.0. 



/i^ enters importantly in the theoretical predictions of dark matter event rates. Thus, if 
yU^ is reduced, the predicted rates generally increase, and if yU^ is increased, they go down. 
From Eq.(4), we see that the non-universal parameters thus can play an important role, as 
for one set of signs (i.e., 5i,3,4 < 0, ^2 > 0) /i^ will be reduced, and the opposite set will 
increase fj,"^. 

3. COMPARISON WITH DAMA DATA 

We compare in this section, the DAMA data with the theoretical expectations for 
SUGRA models. The analysis is carried out using the constraints arising from the radiative 
breaking of the electro- weak symmetry and the accurate method for the computation of the 
relic densitytiJ. 
(i) mSUGRA model 

In this model, all the 6i are set to zero, and the theory is very restricted. Fig. 1 shows 
the experimental bound of the DAMA data for cr^o_p (dotted lines) vs. m^o for the spin 
independent Xi ~ P cross section, where we have combined the 95%CL bound of DAMA 
with the uncertainty in the Milky Way density p^o (0.7 < ^ < 2.3). (The Nal detector is 
sensitive only to the spin independent interaction.) The solid curves are the maximum and 
minimum theoretical cross section cr^o_p for tan/5 < 30 (as one scans over the allowed part of 
the parameter space). In general, the maximum cross section occurs for the maximum values 
of tanjS, so that the upper solid curve occurs for tanjS = 30. The dashed curve corresponds 
to tanj3 = 20 and the dot-dash curve to tanf] = 10. 

The general behavior of the theoretical curves follow from an interplay between the Xi 
early universe annihilation cross section (leading to current relic density) and the Xi - quark 
cross section in the terrestial DM detectoili3. In the relic density analysis, there are two 
regions in the neutralino annihilation cross section. For m^o ~ 50 — 60, annihilation can 
occur rapidly through s-channel Z and h poles, allowing (and sometimes requiring) nio to be 
large (e.g., itiq = (500 — 1000)GeV) so that the lower bound, Vt^oh? > 0.05, not be violated. 



For m^o ~ 50 — GOGeV, the t-channel annihilation through sfermion poles dominates (as one 
now has 2m^o > rrih, Mz) and one requires tjiq to be small (e.g., tjiq ~ ISOGeV) in order 
that the sfermions are sufficiently light so that there is sufficient annihilation to ensure that 
the upper bound Vt^oh^ < 0.3 not be violated. Thus as m^o is increased past 50 GeV, niQ 
decreases. In contrast, the detector x'J'Quark cross section has a major contribution from 
the s-channel squark pole and so increases as rriQ is reduced. In addition, this cross section 
falls off with increasing neutralino mass. 

The above effects can be seen in Fig. 1. The o'^o_„ cross section rises as m^o increases 
and subsequently falls off with increasing m^a. This leads to a maximum theoretical cross 
section at m^a ~ SSGeV, which fortuitously is in the region where the DAMA detector is 
most sensitive. One can therefore accomodate the DAMA data for tanj3 > 8. 
(ii) Models with Non-universal Soft Breaking Non-universal soft breaking can sig- 
nificantly effect the analysis, since as discussed above, /i^ of Eq.(4) can be decreased or 
increased depending on the signs of the 6i. Thus for the cases where 61,62,6^ < and 
64 > one finds that /i^ decreases, and a smaller /i^ tends to enhance the event rates and 
the Xi^P cross-section. We study this case in Fig. 2 where we consider 61 = —1 = —62 and 
S^ = —1 = S4. Here we find that for tanf3 < 30 theoretical predictions often exceed the 
upper limit of the 2cr corridor of DAMA whereas for the mSUGRA case the upper limit of 
theoretical predictions lie within the DAMA corridor (see Fig.l). Similarly the theoretical 
predictions for the case tan(3 < 10 lie significantly higher than the corresponding predic- 
tions for the mSUGRA case. Infact in this case one can achieve consistency with DAMA for 
values of tan/S as low as 6. Further the allowed Xi mass range consistent with the DAMA 
modulation signal extends over a somewhat larger domain of tanf3 relative to the mSUGRA 
case as may be seen by comparison of Fig.l and Fig. 2. 



4. CONSTRAINTS OF PROTON DECAY 

Many SUGRA models which possess a neutrahno dark matter candidate, also give rise 
to proton decay. While predictions of proton decay are more model dependent than in other 
phenomena, it has been observed that cosmological constraints strongly affect predictions 
for the proton lifetime r^, and correspondingly, experimental bounds on Tp affect expected 
DM detector ratesH. 

We consider here GUT groups which contain an SU(5) subgroup with matter embedded 
into the 10 and 5 of SU(5) in the usual way. For these models, proton decay proceeds mainly 
through the mode p — > z/ + K^, with a u and d quark converted into di and ttj squarks (i 
= generation index) by a t-channel chargino (x- ), and the squarks are then converted into 
a z/ and a s-quark by the B and L violating interactions of the superheavy color triplet 
Higgsinos, H^ and i^3ll2lEJ. The proton lifetime is then 

r-i = Tip ^uK)= Y. r(p ^ l^^K^) (5) 

In general, the second generation dominates and in order to get the maximum lifetime, Tmax, 
we assume that the third generation, which enters with arbitrary phase and is about a 20% 
effect, interferes destructively with the second generation. We will also limit the if 3 mass to 
obey Mh-j, < 10 Mq, as larger values of Mh^ would be in the domain of strong gravitational 
effects not being considered in SUGRA models. 

The current Super Kamiokande bound on Tp for the p — >■ PK mode isi^ 

Tp{p -^ uK) > 5.5 X 10^^ yr; 90% C.L. (6) 

One may obtain a qualitative picture of the dependence of Tp on the SUSY parameters 

by noting that the dominant second generation contribution is roughly scaled byOJ 

1 m^o tan/5 „ 
Tip - .A-) ^ ^(^t;^)^ (7) 

Thus large m^o (i.e. large rrig), small rriq and large tan (3 will destabilize the proton. This 
puts strong constraints on the allowed region of the parameter space when combined with 



the relic density constraint 0.05 < Vt^oh? < 0.30. The DAMA data puts additional strong 
constraints on the mSUGRA models. 

In more general situations the Higgs triplet sector of the theory can be more complicated 
and one can have many Higgs triplets which mediate p decay. Thus as discussed in ref.li^ 
in the presence of many Higgs triplets one has a Higgs triplet mass term HiMijHj, and one 
can make a redefinition of fields so that only the Higgs triplet Hi and Hi couple to matter. 
The Higgs triplet couplings to matter then are of the form HiJ + KHi, where J and K 
are bilinear in matter (i.e., quark and lepton fields). If one integrates out the Higgs triplet 
fields, the resulting baryon and lepton number violating dimension five operator is given 
bytJ W4 = —K(M^^)iiJ. Thus for a 2 x 2 matrix M, a suppression by a factor of ~ 3 (e.g., 
by the choice Mu = — 2M22 = M, M12 = M21 = M) in the proton decay amplitude and a 
supression by a factor of ~ 10 in the proton decay rate is easily obtained. In this case the 
proton lifetime can be significantly enhanced. 

We discuss now the numerical analysis of cr^o_p with the inclusion of the proton lifetime 
constraint. It turns out that with the imposition of the proton lifetime limits, one finds 
no points in the parameter space consistent with DAMA data for the minimal SU(5) case. 
However, a non-trivial region of the parameter space does exist for non-minimal SUGRA 
models where one can get consistency simultaneously with both the DAMA data on an- 
nual modulation and the proton lifetime limits. The non-minimalities can arise in various 
forms. We have already discussed non-minimalities in the soft SUSY breaking parame- 
ters where non-minimality implies going beyond mSUGRA to include non-universal soft 
breaking. However, since proton decay involves GUT physics one finds that there are other 
non-minimalities needed to generate acceptable physics at low energy. Thus the minimal 
SU(5) model does not lead to satisfactory quark and lepton mass matrices. One may mod- 
ify it, however, to generate, e.g. the Georgi-Jarlskog texture at M^, leading to reasonably 
correct quark and lepton masses. This requires additional GUT interactions which give rise 
to an increase in r^ of about a factor of 3-5t3. In addition to this, there are other theoretical 
uncertainties in the calculation of Xp, relating to input parameters, e.g. the three quark 
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matrix element of the protonta, and we estimate an additional uncertainty of a factor of 2-3. 
Thus a factor of 10 enhancement is possible for the p lifetime from these factors. Further 
the structure of GUT physics is of course largly unknown. In general as discussd above one 
can conceive of multi Higgs triplets which participate in mediating p decay. As discussed 
above a factor of ten enhancement in the proton lifetime can occur from this source without 
significant fine tuning. Including the enhancements discussed from effects of textures etc. 
one can generate a total enhancement of a factor of ~ 10^. 

An analysis oi Xi—p cross sections with proton decay constraints including the above en- 
hancement factors is given in Fig. 3 for the case of universal soft SUSY breaking parameters. 
One finds that there is a small region of the parameter space where both the DAMA and 
the proton lifetime constraints are satisfied. The allowed x? niass region is close to the peak 
of the DAMA experiment. Next we carry out the analysis for models with non-universal 
soft SUSY breaking. The model which is favorable for both DAMA and proton stability 
is given in Fig.4. It corresponds to the case 6i = —0.5 = —62- The analysis shows that 
with enhancements factors as discussed above one can achieve consistency with DAMA and 
proton lifetime limit over a broad range of Xi mass and within the range given by DAMA. 
This model produces an interesting range of mass spectra for the SUSY particles. Two of 
the supersymmetric particles which are likely candidates for discovery at colliders are the 
Higgs and the lightest chargino. We exhibit the scatter plot of their masses in this model in 
Fig. 5. Their mass ranges are such that they should be accessible at TeV(33). 

5. CONCLUSION 

In this paper we have given an analysis of the annual modulation signal observed by DAMA 
for the direct detection of dark matter within the framework of supergravity modelsli^l. We 
find that SUGRA models with and without non-universalities can easily accomodate the 
range oi Xi ^ P cross-sections needed to explain the observed signal. However, in SUGRA 
models with grand unification, proton lifetime limits impose additional constraints. Thus for 



9 



the case of the minimal SU(5) mSUGRA model it appears not possible to explain the DAMA 
signal and simultaneously satisfy the current proton lifetime lower limits. However, SUGRA 
models with textures and with more complex Higgs triplet structure can allow one to signif- 
icantly enhance the p decay lifetime. Similarly, SUGRA models with non-universalities also 
affect the proton lifetime as well as Xi^P cross-section. It is shown that within the context 
of these non-minimal SUGRA models one can achieve consistency with the signal and a 
satisfaction of the proton lifetime limits. The sparticle mass spectra is also investigated. It 
is shown that the lightest particles in this model are the neutralino Xij the light chargino 
Xi , and the light Higgs h^ which should all be accessible at the TeV(33) and the LHC. 
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Figure Captions 

Figl. Plot of the maximum and minimum oi Xi ^ P cross-section vs neutralino mass for 
mSUGRA for various ranges of tanbeta when other parameters are allowed to vary over the 
ranges mQ^nig < ITeV and Aq is allowed to vary over the range consistent with electro- 
weak symmetry breaking. The various cases correspond to tanf3 < 30 (solid), tan(3 < 20 
(dashed), and tan(3 < 10 (dot-dashed). The minimum dashed and dot-dashed curves overlap 
the minimum solid curve. The dotted curves give the experimental 95%CL range implied 
by the DAMA annual modulation data combined with the uncertainty in ^. 

Fig2. Plot of the maximum and minimum oi Xi~ P cross-section vs neutralino mass for 
the non-universal SUGRA model with 6i = —1 = —62, ^3 = — 1 = 64 for various ranges of 
tanf3 when other parameters are allowed to vary over the ranges mQ,mg < ITeV and Aq is 
allowed to vary over the range consistent with electro-weak symmetry breaking. The various 
cases correspond to tanjS < 30 (solid), tan/S < 10 (dashed), and tanjS < 8 (dot-dashed). 
The minimum dashed and dot-dashed curves overlap the minimum solid curve over most of 
the neutralino mass range. The horizontal curves are as in Fig.l. 
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Fig3. Plot of the maximum and minimum oi Xi ~ P cross-section vs neutralino mass 
for mSUGRA exhibiting the proton hfetime constraint. The maximum and minimum sohd 
curves are for tan(3 < 30 as in Fig.l. The region enclosed by the long dashed curve and 
the minimum solid line is the region allowed under the proton life time constraint with an 
enhancement factor of 10^ as discussed in the text. The region enclosed by the dot-dashed 
lines on the left and the region enclosed by the dot-dashed line and the minimum solid line to 
the right is the region allowed with proton life time constraint with an enhancement factor 
of 20. 

Fig4. Plot of the maximum and minimum of x? "~ P cross-section vs neutralino mass for 
the non-universal SUGRA model with 6i = —0.5 = —62, and ^3 = = ^4 exhibiting the 
proton lifetime constraint. The maximum and minimum solid curves are for the case when 
tan(3 < 30 but without imposition of proton lifetime constraint. The region enclosed by the 
long dashed curve and the minimum solid line is the region allowed by the proton lifetime 
constraint with an enhancement factor of 10^. The region enclosed by the dot-dashed lines 
is the region allowed by the proton life time constraint with an enhancement factor of 20. 
Fig5. The scatter plot of light chargino mass Wi (filled squares) and the light Higgs h^ (open 
circles) for the case 61 = —0.5 = —62 for the set of points that satisfy the DAMA range as 
shown in Fig. 1 and the proton stability constraint corresponding to the area enclosed by the 
dashed curves in Fig4. 
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Fig.l: mSUGRA 

DAMA: 0.7<5<2.3 
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Fig.2: SUGRA 5 ^-l^-e^, 83=-l=54 

DAMA: 0.7<5<2.3 
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Fig.3 
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Fig.4: 5^=-0.5=-52 
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